Localization of an impurity particle on a boson Mott insulator background 
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We investigate the behavior of a single particle hopping on a three dimensional cubic optical lattice 
in the presence of a Mott insulator of bosons in the same lattice. We calculate the critical interaction 
strength between the impurity and background bosons, beyond which there is bound state (polaron) 
formation. We give exact results in the limit of a perfect Mott insulator, where polaron formation 
is equivalent to impurity localization. We calculate the effects of lattice anisotropy, higher impurity 
bands, and fluctuations of the Mott insulator on the localization threshold. We argue that our 
results can be checked experimentally by RF spectroscopy of impurity particles. 

PACS numbers: 03.75.Lm, 71.10.Fd 



I. INTRODUCTION 



Recent experiments using cold atoms in optical lattices have displayed remarkable versatility, creating highly con- 
trollable,isolated, low temperature environments where ideas about many particle quantum mechanics can be tested 
. P, 0, H, 0, [H, [f| . With improvement over the control of system parameters, as well as advancement of novel mea- 
surement techniques such as noise correlations [7| , it seems plausible that a great variety of models will be realized in 
optical lattices. The precision of optical lattice experiments presents another challenge for theory. Measurements of 
' effects that are beyond the simplest descriptions such as mean field theory, are becoming possible (8j. 

One of the new classes of quantum models realized by cold gas experiments is the mixture of different species of 
I ■ atoms [9j. Boson-boson mixtures as well as boson-fermion mixtures have been created. It is now also possible to 
selectively turn on optical lattice potentials for any of the species forming the mixture fiol ] . 

The possibility of experimental realization of atomic gas mixtures stimulated a lot of theoretical interest [111, EH 
E3, EH HE El) E3- F rom the theory point of view, mixtures are appealing as they can display some qualitatively 
new physics stemming from simple ideas in many body theory. There is possibility of pairing due to mediated 
interactions formation of composite particles similar to molecules (l5j . large counterflows of different species or 
even countersuperfluidi ty |T3| . There are also ideas to simulate random potentials using one species as the disorder 
potential for the others [l6| . Two recent experiments about boson-fermion mixtures have shown that the optical lattice 



| experiments are advanced enough so that all these theoretical ideas may now be tested in the laboratory [181 II 



The parameter space for mixtures is very large, with many possible phases and a complicated phase diagram [l5Lll7|. 
The phase diagram of mixtures changes remarkably depending on the strength and type of optical lattices, number 
, of components, nature of interactions, and the identity of the particles. While more general investigations of this 
parameter space, which point out novel phases, are very useful, in this paper, we concentrate on a simple limit and 
, provide some exact results. Experimentally, it would be easy to check such results and gain a better understanding 
CO of the various effects that may be present in the system. 

We investigate the limit where one bosonic species interacts with a single particle of another species. While the 
identity of the external particle does not matter in what follows, we refer to it as the fermion for convenience. All 
the results of the paper are valid for a bosonic impurity as well. We assume that both species share the same 
lattice potential and that the lattice potential is deep enough so that only one band of the lattice is populated. The 
Hubbard-type Hamiltonian [20] for this system can be written as 



^ + IT 

<i,j> i i <i,j> 



H = - f bY, b t b i + -f- E - !) - & E n >- - */ E ftii + U »f E /*. (!) 



where b^,b and f^,f are the bosonic and impurity creation and annihilation operators respectively; rii — b\bi is the 
on-site number operator for the bosons and < i,j > represents a sum over nearest neighbors. The strength of the 
tunnelling terms are characterized by hopping matrix elements tb for bosons and tf for the fermion. Ubb and Ubf are 
the on-site interaction strengths between bosons and between a boson and a fermion respectively. 
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Depending on the ratio between the the tunnelling strength and the interaction for bosons, the system may either 
form a superfluid (SF) or a Mott insulator (MI) state [2l[. While the localization problem can be studied in both 
regimes, we can provide exact results only for the Mott insulator case. In this paper, we concentrate on the problem 
in the Mott regime as the problem is qualitatively different from that of a superfluid interacting with an impurity. 
Thus, we consider tb/Ubb <C 1, give exact results only in the limit tb/Ubb — and calculate corrections due to finite 
boson hopping. For a Mott insulator with no bosons per site, the chemical potential is constrained to 

U b b(n - 1) < Hb < UbbriQ. (2) 

With these considerations for bosons, the fermion can show two qualitatively different behaviors. It can either 
behave like a free fermion with its wavefunction stretching throughout the system, or it may create a defect in the 
Mott insulator and form a bound state with this defect. We calculate the critical interaction strength that separate 
these two regimes. In Section m we study the case of a perfect Mott insulator and obtain the phase diagram for 
localization of the impurity; in Section IIII1 we generalize this exact result to an anisotropic lattice where hopping 
strength in one direction is different from the other two. We consider the effect of higher impurity bands in Section 
IIVI In Section |V1 we calculate the change in the fermion hopping strength due to fluctuations in the Mott insulator 
background. We summarize our results and discuss experimental methods to measure this bound state formation in 
Section ED 



II. LOCALIZATION IN A PERFECT MOTT INSULATOR 



In this section, we calculate the critical interaction strength for bound state (polaron) formation in the limit that 
the boson Mott insulator is perfect, i.e. the hopping strength for bosons is zero (tb/Ubb — 0). We relax this assumption 
in Section [V] and analyze the effects of fluctuations in the Mott background. When boson hopping is neglected the 
Mott insulator background becomes almost inert for the fermion, presenting a spatially independent mean field energy 
shift. In this case the fermion will move with the dispersion relation 



E k =t f 



6 — 2 2_, cos(/cia) 



(3) 



where — 7r/a < ki < ix / a is the crystal momentum of the fermion in the i direction and a is the lattice constant. We 
cannot expect this delocalized behavior of the fermion to continue if the interactions between the background bosons 
and the fermion become very strong. Let us assume that the fermion and the bosons attract each other, i.e. Ubf < 0. 
In this case, it would be energetically favorable to put more bosons at a lattice site and bind the fermion to these 
bosons. This will only happen at a critical interaction strength, beyond which the energy gained by boson-fermion 
attraction is larger than the sum of the kinetic energy cost of localizing the fermion and the interaction energy cost 
of introducing more bosons. 

This critical interaction strength can then be found by investigating the single particle Hamiltonian 

H = -*/ E fUi - Vfofo, (4) 

<ij> 

where we take the site at which the defect is formed as the origin and assume that the defect represents a localized 
attractive potential —V to the fermion. We can then ask for which value of V there will be a bound state in the 
spectrum. This is the discrete version of the problem of existence of a bound state for a localized potential well [12] • 
Just like the continuum version, in one and two dimensional lattices there is a bound state for an infinitesimally small 
attractive potential. In three dimensions, however, there is a certain critical value below which there is no bound 
state. In the case of a finite attractive potential, the dispersion relation takes the following implicit form 

<fe [E + 2 cosiha)] = -Vipo, (5) 

i 

where E = E/tf and V — V/tf are the scaled quantities and (j>k = 'Yln^pi^" 1 is the Fourier transform of the 
particle's wave function. We obtain the relation between the binding energy and the attractive potential by taking 
the Fourier transform of Eq. ([5]) 

i r r r de rr de v de z 1 

(6) 



V }-J-J-, (2tt)3 (6- 2 E< cos + e' 
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FIG. 1: Binding energy e of the impurity as a function of V/tj. The critical interaction strength where the localization begins 
can also be obtained from the figure, i.e. e = for V/tf < 3.96. 

where Oi = fc^a and we take 

E = -6 - e, (7) 

with e > being the binding energy. When e = 0, i.e. at the localization threshold, the above integral can be 
evaluated exactly [23[ and the critical value at which the localization takes place is 

V c = ? = 3.95678, (8) 

(18 + 12V2 - 10V3 - 7V6) [%K(koj\ 

where k$ = [(2 — v / 3)(v / 3 — a/2)] 2 and K is the complete elliptic integral of the first kind. For nonzero e, the integral 
was evaluated by Joyce [24]. Using this result, the potential is obtained as 



where T) = -16z(\/T — 1 + ^/l-9z)- 2 , z = 1/us 2 = 1/(3 + e/2) 2 , and k 2 ± = \ 1 ± ry^/l - \j] - (l - |ry) ^/T^T] 

In the limit of large binding energy, e >> 1, we obtain a linear relation V oc e, which is expected as the particle is 
strongly localized at a single lattice site. 

The exact evaluation of the integral above allows us to calculate not only the critical boundary but also the binding 
energy e(V/tf) of the bound state (Fig. [1]). We now translate our one particle results to the many particle case. Let 
us first assume that the fermion-boson interaction is attractive, Ubf < 0. In this case the simplest defect would be to 
introduce one more boson, thus the attractive potential seen by the fermion will be V = \Ubf\. However, to introduce 
one more boson to a Mott insulator with no particles per lattice site would cost energy, Ubbno — ^ ■ Thus, the phase 
boundary between the free fermion state and the bound state of the fermion and one bosonic defect (polaron) is given 
by 

£ ( Z^£ ) = Ubbn ° ~ lL (10) 

This is by no means the only defect that can be created in the Mott insulator. If the boson-fermion attraction is 
strong enough, it becomes energetically favorable to attract more bosons and form a bound state of two bosons and 
one fermion. The phase boundary for such a defect can be decided by comparing the energy of this state with the 
energy of the bound state of one boson and one fermion. Thus, the equation for phase boundary is 

-U bf ^ _ Ubbnp - \i _ ^ -2U bf ^ _ U b b(2n + 1) - 2/i 
t l t f t f t f 
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FIG. 2: Phase diagram for /i = (no — l/2)Ubb- Numbers in each region show how many extra particles (Ubf < 0) or holes 
(Ubf > 0) are attracted to the localization site. The region marked as > 7 contains all the phases with seven or more extra 
bosons (holes). Phase diagram for this value of /j is symmetric around Ubf — 0. For small boson-boson repulsion Ubb, even for 
small \Ubf\ values, large number of bosons are attracted. While this phase diagram is independent of no, the number of holes 
that are attracted is limited by no. 

One can similarly find the boundaries for bound states with higher number of bosons. 

Another kind of defect would present itself for repulsive interactions, Ubf > 0. For sufficiently strong repulsive 
interactions it would be preferable to create a hole in the Mott insulator state and bind the fermion to this hole. The 
corresponding phase boundary is given by 



Similar to the attractive interactions, it is possible to form bound states of the fermion with more holes. One can 
continue to deplete the Mott state until all the no bosons are removed from the defect site. After this point it would 
be preferable to deplete bosons from the neighboring sites. We have, however, not included such states in our phase 
diagram. In Figs. [2l[3l andHJ we present three phase diagrams for different chemical potentials. Fig. |2] indicates that 
when fi — (hq — l/2)[/ b& , phase diagram is symmetric around Ubf = 0. This is expected as this value of ji corresponds 
to lobe centers of the Bose-Hubbard phase diagram where there is particle-hole symmetry. One can also notice from 
this diagram that when Ubb is close to zero, even for small \Ubf\ values, the fermion can be bound to a large number 
of bosons. If we take \i = (n — 1/4) [4b as in Fig. [3l the symmetry around Ubf = is broken and for the repulsive 
interactions it is harder to attract holes. Fig. [4] represents the opposite case (/i = (no — 3/4) [/(,(,) where stronger 
interactions are required to attract particles. 

We believe that the phase diagram can be checked experimentally. While it would be possible to modify Ubf by 
an interspecies Feshbach resonance, an easier route would be to change tf which is controlled by the strength of the 
optical lattice. The localized impurity states can be distinguished from free fermion states as their mean field shifts 
would be different; in principle, RF spectroscopy [H|,[2(| would directly detect the difference in the mean field shift. 
Although the calculation was carried out for a single impurity, we expect these results to be quantitatively correct for 
a small density of fcrmionic impurities over a bosonic Mott insulator background. Essentially, if the inverse of Fermi 
momentum is much larger than the lattice spacing, then the fermions would hardly effect each other's behavior. 

In this section, we obtained the phase diagram for the interaction of one impurity particle with a perfect Mott 
insulator. In the next three sections, we investigate how this ideal situation is affected by lattice anisotropy, higher 
impurity bands, and fluctuations of the Mott insulator. 




Ubbjno - 1) + /z 



(12) 



III. EFFECTS OF LATTICE ANISOTROPY 



In the optical lattice experiments, it is possible to change the strength of the laser beams forming the lattice, 
hence realize a model system where the lattice is not isotropic. For a quantitative comparison with experiment, it is 
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FIG. 3: Phase diagram for a — (no — l/4)Ubb- Symmetry in Fig. [2] is broken and particle attraction is easier than the hole 
attraction, since the chemical potential is increased with respect to the symmetry point. To attract a hole one needs higher 
boson-fermion interaction for the same Ubb- 




FIG. 4: Phase diagram for /i = (no — 3/4)Ubb- As opposed to Fig. [3] to attract a particle one needs higher boson-fermion 
interaction. 

necessary to take this effect into account. We assume that the hopping strength for the fermion is different in one 
direction compared to the other two directions and calculate the effect of such anisotropy on the phase diagram of 
the previous section. 

The localization threshold for the anisotropic case can also be calculated analytically. Thus, in the following 
discussion we need not assume that the anisotropy of the lattice is small. As a simple limit, we obtain the two 
dimensional lattice localization problem as the layers are decoupled. For that limit, even the smallest attraction 
between the fermion and bosons causes bound state formation. As the coupling between the two dimensional layers 
is increased the critical interaction needed for bound state formation monotonically increases. 

Because of the anisotropy, the single particle Hamiltonian in Eq. is modified as 

H = -*/ E fifi - */ E flff - Vf o^ ( 13 ) 

<i,j> <i,j> 
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where we take the hopping term in the z direction to be t'* 7^ t/. We obtain the relation between V and e as 



1 

V 



7T P 7T p7T 



7T •/ — 7T «/ — 7T 



d8 x d6 v d6 z 



(2tt)3 [4 - 2 £ l=a; cos 0i + 2r(l - cos 6 Z )] + e' 



(14) 



where E = —4 — 2t — e and r = t'^/tf. For e = this integral can be evaluated exactly [23|. The critical value for the 
localization is found to be 



^(V2VTTT- V2+T) (f ) 2 K[fc + (r)]^[fc_(r)] ' 



where 



Mr) 2 = 



- ( V2 VTT7 - V2T7) ( V2T7 ± V2) 



(15) 



(16) 



As t increases, i.e. the anisotropy of the lattice increases, the critical value for the potential increases and the 
localization becomes more difficult (Fig. [SJ. As tj — > 0, the system becomes two dimensional and there is no 



threshold for localization, as expected. For large t'f/tf, V c ~ «/iyt/, which gives V c — » in the one dimensional limit, 
tf — ► 0. Moreover, it is possible to evaluate the integral for nonzero e [HI, yielding 



where w — 2 + r + e/2, z = 1/u; 2 and 



1 1 



2(V1-(2-t)2z + v /1-(2 + t)2z) 
W (|) 2 ^[A: + (r)]^[fc_(r)] 



(17) 



£^ — - - - 

k± ~ 2 2 



v /1-(2-t) 2 z + v /1-(2 + t) 2 z 



(18) 



y/l + (2- t) Vly/l - (2 + r)VI + 1 - (2 - r)v^/ 1 + (2 + r)v/z 



±16z+ VI - r 2 ; 



1 + (2 - TlVW 1 + ( 2 + + V 1 - ( 2 - ^VzV 1 - ( 2 + r 



Using this exact result, the phase diagram can be obtained for arbitrary t. In Fig. [6l we display the phase diagram 
for r = 1.5. Comparing Fig. [S]with Fig. [2] (isotropic case) we see that the phase boundaries are closer to the Ubf 
axis. 



IV. EFFECTS OF HIGHER IMPURITY BANDS 



An important point one always has to keep in mind about the optical lattice experiments is that the effective 
Hubbard models, such as Eq. Q, are obtained by projecting the system into the lowest band of the lattice [2fJ. This 
procedure is expected to describe the low energy physics as long as the band gaps are larger than the temperature 
and interaction scales in the problem. In the equivalent language of Wannier functions, this condition corresponds to 
requiring the Wannier function of each lattice site to be undisturbed by interactions. 

In the context of the current problem, we discussed the critical hopping strength that is needed to localize the 
impurity particle to a small region, which is of the order of one lattice site. The precise determination of the Hubbard 
model parameters such as Uf,f depends on the microscopic model one starts from. For the Hubbard model to work 
correctly, the Wannier functions for the impurity must be unchanged even if the impurity particle is localized to one 
lattice site. As a localized impurity attracts (or repels) extra particles (holes) to its localization site, one may expect 
the on-site wave function of the localized particle to be different from the Wannier functions at other lattice sites. 
This is essentially considering the coupling of the localized particle to higher impurity bands, and should be a small 
effect controlled by the parameter where A/ is the width of the first band gap of the impurity bands. Thus, the 
effect we are considering in this section would be important only if the impurity particle is highly mobile in the lattice, 
while the interaction between the background particles and the impurity is strong enough to localize the particle (Fig. 

ED. 
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FIG. 5: The critical value of interaction V c /tf as a function of lattice anisotropy characterized by t'f/tf (Eq. H5J)). If the 
hopping parameter t'f increases, anisotropy of the lattice increases and the localization becomes more difficult, r = t'f/tf = 1 
gives V c /tf for the isotropic case. 




u bb /t f 

FIG. 6: Phase diagram in the presence of lattice anisotropy (for r = t'f/tf — 1.5 and /i = (no — l/2)Utb)- To be compared with 
Fig. [2](r = 1, fi = (no — 1/2) [/(,(,). One can see that anisotropy with r > 1 causes the localization threshold to move to higher 
values of |C4/|. 




FIG. 7: Schematic representation of the effect of higher impurity bands to the hopping parameter. If the localized impurity 
attracts extra particles (holes) to the localization site, the local wave function of the impurity particle changes. Then the 
hopping parameter for this site is different from that for the other sites. 
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FIG. 8: The critical value of interaction V c /tf as a function of r = tf/tf (Eq. (|22p ). As the ratio of the hopping parameters r 
increases, localization occurs for smaller values of the interaction. 

In such a case, the system can still be modelled by a Hubbard model where the hopping strength between the 
localization site and its neighbors (t'j) is different from the hopping strength between any other neighboring sites in 
the lattice (tf). These hopping strengths can once again be calculated by looking at the overlaps of the localized 
wavefunctions between neighboring lattice sites (2pj |. 

In this case, we take the single particle Hamiltonian as 

H = J2 fUi - (*/ - */) E flfrn(Sl + S m0 ) - Vftfo. (19) 

<i,j> <i,m> 

Calculations similar to those performed in the previous sections yield 

i _ r* pr de x d0 y de z 2(t-1)^.cos0; 

i-ixi-ixJ-ix (2ir) 3 6-2 ^J. cos 6»i+e 
^ = pr pr /-7T d6xd9 v d0 z 1 ' (^0) 

J -ix J -ix J -t; (27r) 3 6-2^ .cos6' i +c 

where -B = — 6 — e. Evaluating this integral exactly [27|, we obtain the relation 



where = -16z( A /T _ I + VI - 9z)- 2 , z = = 1/(3 + e/2) 2 , and k 2 ± = \ 1 ± 77WI - \i) - (l - |ry) VT"^ 

When e = 0, the critical value for the potential is found to be 

V c » 3.95678[1 - 0.51622(t - 1)]. (22) 

As can be seen in Fig. [51 the coupling to higher impurity bands can substantially change the critical value for 
localization. If Ubf < 0, we expect a narrowing of the local wave function (as in Fig. [7]), then t'^ < tf and consequently 
localization is harder V c (r) > V c (t = 1). Similarly if Ubf > we expect easier localization. In general one would then 
expect each different polaron state to have a different r value. Still to gain a basic understanding of this effect we 
obtain the phase diagram (Fig. [9|) using constant value of r = 1.5. 

V. EFFECTS OF FLUCTUATIONS OF THE MOTT INSULATOR 



The ideal Mott insulator state is achieved only when the boson hopping term in the Bose-Hubbard Hamiltonian 
is zero. When there is a small but nonzero hopping probability for the bosons, the ground state of the MI contains 
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FIG. 9: Phase diagram obtained when the effect of higher impurity bands is taken into account (r = 1.5, [i = (no — l/2)f&f>). 
This effect is modelled by the parameter t, which is the ratio of the hopping strength between the localization site and its 
neighbors to the one between any other neighboring sites. Compare this figure with Fig. [2] (r = 1, n — (no — l/2)Ubt)- One 
can see that localization is easier if r > 1. 



virtual particle and hole excitations. When tb/Ubb gets larger, these excitations gain more amplitude and will finally 
destroy the order in the MI, and cause a transition to the SF state. 

In the previous sections, we considered the MI to be a perfect insulator by setting boson hopping to zero. In this 
section, we consider the effects of small but nonzero hopping on our calculations. We always work in the limit that 
the boson hopping th is the smallest energy scale of the system and require that the system be away from the MI-SF 
transition boundaries. Under these conditions, we can use perturbation theory to investigate the effects of particle 
and hole fluctuations on the localization of the impurity particle. 

There are two main consequences of turning on the boson hopping. First, the compound object formed by the 
impurity and the extra bosons (polaron) will become mobile. The width of the polaron band will be proportional to 
tft^, where S is the number of extra bosons (or holes) forming the polaron. As th is the smallest energy scale in the 
problem, this polaron mobility effect will be small, especially for polarons containing more than one boson or hole. 
In other words as tf ^> %, the fermion is much more mobile than the boson (or the combined polaron). The mobility 
of the polaron will then have a small effect on the localization problem considered above. 

The second effect is the change of effective mass (or effective hopping tf) of the impurity particle due to scattering 
from the fluctuations of the MI. The first contribution to such processes will be proportional to (Ubf/Uhb) 2 for the 
fermion to scatter from a particle-hole pair, and (tb/Ubb) 2 for the pair to be excited. Because of the presence of U 2 ^ 
which is large in the limit that we are interested in, this second effect will be dominant over polaron mobility. We 
should also mention that it is only lowest orders of the perturbation theory that these two effects can be considered 
separate processes. 

To the lowest order in tb/Ubb the boson wavefunction can be written as: | \f r ) =| MI) + X)<j j> tybj \ MI). The 
state of the system without boson- fermion interaction is then | ty) | k) , where | k) is the wavefunction of a fermion 
moving with lattice momentum k. 

By treating the boson-fermion interaction as a perturbation 

V = U if Y,mftfu (23) 

i 

we can write the second order energy shift as 



A(2) _ ^ | (k' | (MI | bjb p V[ \MI) + ^ E <tj> b\ bj | MI)] | k) | 2 

~Jfc e h -<p» + e„) > (24) 



where the only non- vanishing term is 
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(k' I (MI | bl b V^- J2 bib, | MI) | k) = ^^fno(n + 1) [e«*-*>^ _ e <C*-S>-«.]. (25 ) 

<Ab ' <^bb 

<i,3> 

Again considering the lowest non-vanishing order contribution we get: 

U bb N <a y >ki £k -(u bb + £k ,) 

- & ± n i (m + f d3fc/ 12-4E,cos[(^-fc) a ] 



'bb 
^bb 



12-28-^ Vcos^) , ^ L «l. (28) 



C/bb 



Hence, we see that the first effect of the fluctuations of the MI is to renormalize the hopping strength of the impurity 
particle to 

*/^/(l-14^K + l) 2 ^§^. (29) 

Essentially a decreasing tf means the effective mass of the particle is larger, and it becomes easier to localize the 
particle. Thus, the effects of boson hopping on all the situations discussed in previous sections can be obtained by 
scaling the phase diagrams by the renormalized value of tf . 

VI. CONCLUSIONS 

We study the localization of a single impurity particle in a lattice containing a Mott insulator of bosons. This is a 
simple limit of the complex physics of the Bosc-Hubbard model for mixtures of different atomic species. The impurity 
particle has two distinct types of behavior; it may either move freely throughout the lattice or may choose to localize 
at a certain lattice site by attracting extra bosons or holes. In the limit of a perfect Mott insulator, we calculate the 
boundary between these two phases as well as the number of extra bosons (holes) forming the bound state exactly. 
Our result for the phase diagram is given in Fig. [21 

We believe that this phase diagram can be checked experimentally. In recent experiments boson-fermion mixtures 
were created in the parameter regimes that we consider in this paper (lR [l9| . If the density of fermions is small 
enough, one can disregard the many-particle nature of the fermions and consider the localization problem for one of 
them. To determine whether a fermion is localized, RF spectroscopy would be an ideal tool. The mean-field shift of 
a localized fermion would be larger as it sees more background bosons on the average. Equivalcntly, one can say that 
the binding energy of the fermion calculated in Section |TT] would be directly reflected in the RF spectra for fermions. 
Different polaron states, containing multiple numbers of bosons or holes can be distinguished similarly. 

The calculations presented in this paper were carried out with the assumption of a homogenous infinite system. In 
cold gas experiments there is always a confining trap, which in general complicates the correspondence between the 
infinite system predictions and experimental results. However, when a lattice boson system is driven deep into the MI 
regime, the density profile of the system consists of MI plateaus separated by thin SF regions. Within each MI region, 
the fermion would see a flat interaction potential, and when the fermion is localized, the width of the wavefunction 
of the fermion becomes of the order of few lattice sites. Thus, we expect our calculations to accurately reflect the 
transition boundaries, especially for transitions between different polaron states. 

After the ideal case, we consider three effects which may play a role in the experiments. The first case we consider 
is the possibility of tunnelling anisotropy. As the strength of the lattice is determined by the laser intensity (and 
beam waist), experiments occasionally have such anisotropy. We generalize the exact results to this case and find 
that it becomes easier to localize the impurity when the system becomes more two dimensional, as expected. Next, 
we consider the effect of higher bands of the impurity particle. We argue that this effect can be taken into account 
by modifying the hopping strengths between the localization site and its neighbors, and obtain the phase diagram. 
Finally, we consider the effect of small but finite tunnelling strength for the bosons forming the MI, and calculate the 
effective hopping strength for the impurity particle in the presence of fluctuations. 

We believe that our exact results about impurity localization on a Mott insulator background provide a starting 
point for the investigation of the complex phase diagram of mixtures in optical lattices. 
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